'''
Created on 21/04/2012

@author: nicolas
'''

import pylab as p

dTime = 0.001           #intervalo de muestreo
NFFT = 512              #largo de ventana
Fs = int(1.0 / dTime)   #frecuencia de muestreo
noverlap = 16

#defino el tiempo de muestreo
t = p.arange(0, 1, dTime)

#defino las funciones
def sin1(t):
    return p.sin(2 * p.pi * 100 * t)

def sin2(t):
    return p.sin(2 * p.pi * 150 * t)

#defino la suma de senos
def sinSum(t):
    return (2 * sin1(t) + sin2(t))

#realizo el muestreo
x = []
x2 = []
x3 = []

for i in t:
    x.append(sin1(i))
    x2.append(sin2(i))
    x3.append(sinSum(i))

# Grafico
p.figure(1)

p.subplot(211)
p.title('100 hz')
p.plot(t, x)

p.subplot(212)
p.title('150 hz')
p.plot(t, x2)

p.savefig('punto_1_sins.png')

p.figure(2)
p.title('Espectograma')
p.specgram(x3, NFFT=NFFT, Fs=Fs, noverlap=noverlap)

p.savefig('punto_1_spectogram.png')

p.show()
